a-clustering effects in dissipative 12 C+ 12 C reactions at 95 MeV 
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Dissipative 12 C+ 12 C reactions at 95 MeV are fully detected in charge with the GARFIELD and 
RCo apparatuses at LNL. A comparison to a dedicated Hauser-Feshbach calculation allows to select 
events which correspond, to a large extent, to the statistical evaporation of highly excited Mg, as 
well as to extract information on the isotopic distribution of the evaporation residues in coincidence 
with their complete evaporation chain. Residual deviations from a statistical behaviour are observed 
in a yields and attributed to the persistence of cluster correlations well above the 24 Mg threshold 
for 6 a's decay. 
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Since the first heuristic proposition of a-chains as pos- 
sible building blocks of even-even nuclei in the late six- 
ties the subject of a-clustering has been a central 
issue in nuclear physics and has even witnessed a gain 
of interest in recent years On the theoretical side, 
highly sophisticated ab-initio calculations have shown 
pronounced cluster features in the ground state of a large 
number of light nuclei Q as well as in some excited states 
around the threshold energy of breakup into constituent 
clusters, showing that cluster correlations are indeed an 
ubiquitous feature of quantum few-body systems down 
to the femtometer scale. Concerning experimental re- 
search, rotational bands consistent with a-cluster struc- 
tures have been identified in different even-even light 
nuclei and shown to persist even along their isotopic 
chains 0]. Exotic non-statistical decays of these cor- 
related states have been evidenced in the recent litera- 
ture 0]. 

A natural extension of the concept of nuclear clusters 
concerns nuclear molecules. Molecular states have been 
seeked for in nuclear reactions since the early days of 
heavy-ion science. In particular, several interesting res- 
onances have been observed in the 12 C+ 12 C reaction in 
the inelastic Q and a-transfer channels @ . These studies 
suggest that resonant structures persist in the 24 Mg sys- 
tem up to around 50 MeV excitation energy, a remarkable 
result as a pure statistical behaviour might be expected 
due to the extremely high number of available states at 
such high excitation. Concerning the a-transfer channel, 
experimental results have been reproduced by coupled 
cluster calculations 7] where the cross section is domi- 
nated by a four-cluster (a + a) + (a+ 12 C) state of highly 



excited 24 Mg around 30 MeV. Because of the remarkable 
persistence of cluster correlations at high excitation ener- 
gies, the question naturally arises whether such correla- 
tions might affect other more dissipative channels which 
are typically associated to the formation of a compound 
nucleus, that is a system whose decay is assumed to be 
fully decoupled from the reaction entrance channel and 
governed by purely statistical laws. 

This effect might be experimentally seen as an excess 
of cluster production with respect to the prediction of 
the statistical model, provided that the ingredients of 
the latter are sufficiently constrained via experimental 
data. To this aim, we have performed an exclusive and 
(quasi) complete detection of the different decay prod- 
ucts emitted in 12 C+ 12 C dissipative reactions at 95 MeV, 
and compared experimental data to a dedicated Hauser- 
Feshbach code, with transmission coefficients and level 
densities optimized on the A w 20 region [8| . 

We will show that all the observables of dissipative 
events are fully compatible with a standard statistical 
behaviour, with the exception of a-yields in coincidence 
with oxygen residues. The observed anomalies are ten- 
tatively attributed to clustering effects which appear to 
survive even in the most dissipative events. 

The experiment was performed at the LNL (Labora- 
tori Nazionali di Legnaro), with the 12 C beam provided 
by the XTU TANDEM accelerator. The experimental 
setup was composed by the GARFIELD detector, cov- 
ering almost completely the polar angular range from 
30° to 150°, and the Ring-Counter (RCo) annular de- 
tector, centered at 0° with respect to the beam direc- 
tion and covering forward laboratory angles in the range 
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5° < < 18°. Details on the apparatus and the experi- 
mental techniques can be found in Ref. 0, [l(| • The com- 
bination of the two devices allows a nearly-47r coverage 
of the solid angle, which, combined with a high granular- 
ity, permits to measure the charge, the energy and the 
emission angles of nearly all charged reaction products, 
allowing an excellent discrimination of different reaction 
mechanisms. 

Events with a single Z > 2 fragment detected in the 
RCo (identified only in charge) in coincidence with one or 
more light charged particles (LCP) Z < 3 in GARFIELD 
are selected. LCP isotopic identification has been per- 
formed through Fast and Slow Csl analysis [11], with 
energy thresholds 12| similar to other experimental 4tt 
devices employing the same identification method |13j . 
Only events where the total charge of the entrance chan- 
nel Zdet = 12 is detected are kept for the analysis, unless 
specified in the text. We expect that most of the selected 
events should correspond to fusion-evaporation, with the 
residue detected at forward angles in coincidence with 
evaporated particles detected at central center-of-mass 
angles covered by GARFIELD. Experimental data are 
therefore compared to the predictions of a Monte Carlo 
Hauser-Feshbach code [1] for the evaporation of the com- 
pound nucleus 24 Mg, at Ej us — 61.4 MeV, correspond- 
ing to complete fusion source, and filtered through a 
software replica of the experimental set-up. The max- 
imum angular momentum for the fused system is as- 
sumed J n m ax = 12 H, coming from the systematics in 
PACE4 [14J. The comparison of various experimental 
observables and code calculations is used to validate the 
parameterizations of statistical model ingredients imple- 
mented in the code. In particular, an exclusive study 
of completely reconstructed decay chains has been per- 
formed, in order to get an insight into the deviations from 
a statistical behaviour observed in the decay. 

Fig. [1] displays the energy spectra of protons and a 
particles detected in GARFIELD in coincidence with the 
most abundant residues. If not explicitely stated, all dis- 
tributions are shown normalized to a unitary area, so that 
a comparison among their shapes is easily done, and all 
energies are given in the laboratory reference frame. Ex- 
perimental data are always shown with statistical error 
bars, when visible. We can see that a very good reproduc- 
tion of the proton energy spectra is achieved in all chan- 
nels, while a large discrepancy in the shape of the distri- 
bution appears for a particles in coincidence with oxy- 
gen residues. The shape of particle spectra produced by 
a statistical process is determined by the interplay of all 
physical ingredients entering in the decay, notably trans- 
mission coefficients, angular momentum distribution and 
level density parameters. Nevertheless, these different in- 
gredients can be largely disentangled f]~5j | : transmission 
coefficients define the shape of evaporated spectra in the 
Coulomb barrier region; the level density mostly affects 
the slope of the exponential tail; and the angular momen- 
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FIG. 1: (Color online) Proton (upper part) and a (lower part) 
energy spectra in complete Zdet = 12 events detected in coin- 
cidence with a residue of charge Z res , indicated in each figure 
column. Experimental data (dots) are compared to model cal- 
culations (red lines). Spectra are shown for the most abuntant 
residues, 7 < Z res < 10. 



turn steepens the tail of the spectrum of heavy particles 
(a's) with respect to light ones (protons). On the basis of 
such considerations it is found that no common choice on 
the parameters can be done in our calculations in order to 
reproduce at the same time protons and a energy spectra. 
The best model reproduction is obtained with standard 
fiducial values for the transmission coefficients, level den- 
sity parameters and angular momentum distribution Q , 
as shown in Fig. [TJ The very good reproduction achieved 
for protons suggests that the discrepancy found for a par- 
ticles reflects an out-of-equilibrium emission and mostly 
concerns well-defined channels where oxygen is produced. 

To understand the deviation, a more exclusive channel 
by channel comparison is therefore needed. Fig. [5] dis- 
plays the experimental (left panel) and theoretical (right) 
correlation between the energy of the oxygen fragment 
E res and the sum of the energies of the two a-particles 
Ea.i in 12 C( 12 C," 4 0)ao; events. The lines represent the 
kinematical locus Q kin = E res + J2i=i Ea i ~ E beam = 
— 15.78 MeV, where E\, eam is the beam energy. This locus 
divides the {E res ,z2i=i Eai) plane into two regions: in 
the one above, we find two correlated bands centered re- 
spectively at Q kin = (1.27±1.40) MeV and (-5.35±1.70) 
MeV, while a broader region extends below the locus 
up to a high amount of missing energy. In the follow- 
ing, we will refer to these regions as respectively the 
non-dissipative and dissipative events region, and we will 
adopt for them the notation Q > and Q < . In the sta- 
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FIG. 2: (Color online) Experimental (a) and theoretical (b) 
correlation between the energy of the oxygen fragment and the 
sum of the energies of the two a particles in 12 C( 12 C, A 0)aa 
events. The line represent the kinematical locus Qkin = 
Eres + Y^ii=i EoLi — Ebeam = —15.78 MeV at the opening of 
the 4-body channel 12 C( 12 C, 15 0)n«a. 



tistical model interpretation, the two correlated bands 
correspond to a-decay chains, starting from the 24 Mg* 
compound nucleus and leaving an 16 O residue either in 
its ground state {Qkin = —0.11 MeV) or in one of its 
excited bound states which are not resolved in the ex- 
periment, with an average energy expense Qkin ~ —6.55 
MeV. The broadening observed for the experimental Qkin 
bands reflects the energy resolution of our setup for the 
energy balance in triple coincidence measurements, which 
adds up to AE/E < 4%. The locus Q Hn = -15.78 MeV 
gives the threshold Q-value where the 4-body channel 
12 C( 12 C, 15 0)rtaa opens. In terms of energy, this is the 
less expensive channel involving an oxygen residue and 
particles other than a's. Neutrons are not detected in 
this experiment, and the broader distribution observed 
for lower Qkin values is due to events in which neu- 
tron(s) emission has taken place, and their kinetic energy 
has not been collected. According to model calculations, 
the (n,2a, 15 0) channel should absorb the largest cross- 
section in this region. 

Besides the common pattern observed for the (2a, A 0) 
energy correlation, a difference in the relative population 
of the regions is evident between experimental data 
and calculations in Fig. [2] In particular, a much higher 
percentage of (2a, A 0) events populates the Q > region 
in the experimental sample: (37±5)%, with respect to 
(9±1)% according to model predictions. Also in the 
following, when giving experimental percentages, the 
associated error is an upper limit taking into account 
both the statistical error and the possible 3 He-a con- 
tamination. Events falling in the Q > region correspond 
to exit channels with low energy dissipation. Therefore, 
this region can also be populated by reactions not 
proceeding through an intermediate compound nucleus 
state. The larger experimental branching ratio for the 
(2a, 16 0) exit channel in the Q> region probably reflects 
a contamination of direct (a-transfer/pick-up) reactions, 



in competion with fusion-evaporation. 

Fig. [3]a displays the comparison between data and cal- 
culations for the energy spectrum of a particles detected 
in coincidence with an oxygen residue and emitted in dis- 
sipative events (full dots). As a reference the a energy 
distribution without any Qkin selection (empty dots, see 
Fig. [T]) is also plotted. A great improvement in the agree- 
ment between data and calculation is achieved when we 
limit the comparison to the Q< event class. Nevertheless, 
a residual discrepancy is observed, and the experimental 
distribution is still not reproduced in its shape. To in- 
vestigate the origin of this deviation, in the middle (6) 
and right (c) panels of Fig. [3] we show the a energy 
distributions associated to channels where the oxygen 
residue is in coincidence with, respectively, two or only 
one a('s). We can see that the shape of the energy spec- 
tra depends on the channel. The statistical calculation is 
able to satisfactorily reproduce separately the energy dis- 
tributions of a particles for dissipative 2a's and la com- 
pletely reconstructed decay chains, while the sum of the 
two processes is not correctly reproduced. Experimen- 
tally, (63±5)% of the total cross section for completely 
reconstructed dissipative decays with a a particle and an 
oxygen fragment in the outgoing channel is absorbed by 
(2a, 16 0*) channels, where the oxygen is excited above its 
neutron emission threshold, while according to the the- 
oretical predictions these channels should represent only 
(10±1)% of this class of events. In conclusion, the resid- 
ual deviation observed in the leftmost panel of Fig. [3] ap- 
pears to originate from a branching ratio discrepancy for 
the (n,2a, 15 0) (and, with a minor contribution, for the 
2n,2a, 14 0)) decay channel with respect to decay chan- 
nels involving a single a emission, even if a contribution 
from high energy 7 emission from collective states cannot 
be a priori excluded. 

This finding indicates that cluster correlations associ- 
ated to the 12 C+ 12 C system persist up to higher center- 
of-mass energies than previously expected from the anal- 



ysis of the inelastic channel [16[, and lead to a non- 
statistical behaviour in the decay of the highly excited 
24 Mg. Given the high excitation energy E* = 61.4 MeV 
it is not possible to associate the extra yield to a sin- 
gle isolated state with a well defined angular momentum 
and parity. Nevertheless, interesting information can be 
extracted from Fig. [H which displays the a — a relative 
energy distribution for the Q< and the Q > region, in 
comparison with the statistical model calculation. We 
can see that a clear peak appears in the most dissipative 
events, which can be associated to doorway ( 8 Be- 16 0*) 
and ( 9 Be*- 15 0) states, where the excitation energy of the 
16 ( 9 Be) is above the neutron emission threshold, and 
which is not present in the statistical model. The ab- 
sence of a 8 Be resonance in the less dissipative events 
suggests that these reactions do not originate from a 
(a + a)+(a+ 12 C) doorway state as it was suggested to 
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FIG. 3: (Color online) For dissipative events (as defined in 
the text) with Zdet = 12, experimental (black dots) and cal- 
culated (red lines) a energy spectra in coincidence with an 
oxygen are compared for: all channels involving at least la 
(a) and, separately, channels involving 2a's (b) and only 1 a 
(c). As a reference, in the leftmost panel, also the a energy 
distribution without any Qkin selection (same as in Fig. [TJ is 
shown (empty dots). 
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FIG. 4: (Color online) Experimental (black dots) and calcu- 
lated (red lines) relative energy distributions of the two a's 
found in coincidence with an oxygen in dissipative (left) and 
non-dissipative (right) events. In the left panel, a zoom on 
the low relative energy region is shown in the figure inset for 
experimental data. 



be the case at lower energies Q. 

Since we have interpreted the observed extra ct-yield 
from dissipative events as associated to the production of 
neutron-poor oxygen isotopes, we expect that a signature 
of the anomalous branching ratio should be evidenced 
in the isotopic residue distribution. Unfortunately the 
residue mass is not directly accessible in this experiment, 
because of the low energy of heavier fragments. However, 
due to the completeness of the detection, this informa- 
tion can be approximately deduced from the calorimetric 
energy distribution, as explained hereafter. The excita- 
tion energy Ej us can be estimated event- by- event from 
the calorimetric energy balance E* ,(A res ), provided an 



hypothesis is made on the residue mass number A res : 

Nc 

E* cal {A res ) = £ Ef M + N n (Ares) ■ (E% M ) + Q(A res ) 

i 

(i) 

where N c ( N n ) and Ef M ((E^ M )) are respectively 
the charged products (neutron) multiplicity and their 
center-of-mass kinetic energies, and Q is the decay Q- 
value. Indeed, the emitted neutron number is uniquely 
determined by mass conservation and their energy can 
be (in average) estimated from the measured proton en- 
ergy with the subtraction of the Coulomb barrier. An 
estimation of the mass of the residue can thus be ob- 
tained by minimizing in each event completely detected 
in charge the quantity Se* al (A res ) = \E* al (A res ) - Ef us \. 
Fig. [SJa displays the resulting estimated oxygen isotopic 
distribution. As expected, we find in the data an ex- 
tra production of neutron-poor isotopes, coherent with 
the finding of an extra experimental cross section for the 
(xn,2a, 16 ~ x O) channel. The reconstructed experimental 
isotopic distribution for the considered residue is found 
to be generically wider with respect to model predictions. 

Finally, it is important to stress that, with the ex- 
ception of the anomalous branching ratio towards this 
channel, the detected events where the dissipated energy 
overcomes the four-body threshold associated to 15 O for- 
mation (Qkin = —15.78 MeV) are fully compatible with 
a complete fusion pattern followed by the compound nu- 
cleus decay of 24 Mg*. Indeed the Hauser-Feshbach cal- 
culation satisfactorily reproduces all the different observ- 
ables that can be constructed with the exclusion of Q> 
events from the analysis, and that do not directly involve 
a-oxygen coincidences. Fig. [5] displays some selected ex- 
amples, as the velocity distributions of the heaviest frag- 
ment and the charge distribution of reaction products, 
this latter both for complete (Zdet — 12) and almost- 
complete (Zd et > 10) events, in order to show that no 
important bias is induced by the stringent completeness 
conditions. The global quality of the agreement between 
data and calculations makes us claim that we have been 
able to select in our data-set a component which is as 
close as possible to the statistical decay of a hot equi- 
librated source, corresponding to the complete fusion 
24 Mg* source. 

In conclusion, in this paper we have presented an 
exclusive analysis of dissipative 12 C+ 12 C collisions at 
95 MeV, fully reconstructed in charge. A detailed 
comparison to a dedicated Hauser-Feshbach calculation 
shows an abnormally high branching ratio towards 
the (2a, 16 0*) channel with respect to the statistical 
expectation, where the oxygen is excited above its 
neutron emission threshold, which corresponds in part 
to the population of a doorway ( 8 Be- 16 0*) configu- 
ration. This extra yield could be due to a failure of 
the Hauser-Feshbach theory for the decay of highly 
excited 24 Mg, possibly due to the a-correlations of this 
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FIG. 5: (Color online) Panel (a): the experimental oxygen iso- 
topic distribution (reconstructed via the minimization of the 
quantity 5e* cal defined in the text, black dots) is compared to 
the model predictions (red histograms) for all events. Panel 
(b) and lower part, respectively: charge distribution of re- 
action products and velocity distributions of most abundant 
evaporation residues (charge Z res indicated in parenthesis) 
are compared for data (dots) and model calculations (lines) 
after the exclusion of non-dissipative events. For the charge 
distribution, experimental data for complete (Zdet = 12, black 
points) and quasi-complete (Z det > 10, empty points) events 
are compared to model calculations (solid/dashed lines re- 
spectively). Mass and charge distributions are normalized to 
the number of events retained in the analysis. 

nucleus which are already suspected to lead to specific 
high- lying ( 16 0-2a) resonances [17| . Alternatively, the 
compound nucleus hypothesis of full decoupling between 
entrance and exit channel could be questioned due to 
a correlations in the 12 C+ 12 C entrance channel. These 
two hypotheses cannot be discriminated by means of 
the present experimental information, but new data on 
14 N+ 10 B at 80 MeV are currently under analysis to give 
a definitive answer to this issue. 
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